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the simulation of an application of firefly algorithm (FA) for optimally 
locating the most suitable placement and capacity of distributed generation 
Keywords: (DG) in IEEE 33-bus radial distribution network. This strategy aims at 
minimizing losses together with improving the voltage profile in distribution 


Distributed generation network. The losses in real power and voltages at each bus are obtained using 


Firefly algorithm et load flow analysis which was performed on an IEEE 33-bus radial 
IEEE 33-bus distribution distribution network using forward sweep method. The proposed method 
system comprises of simulation of the test system with DG as well as in the absence 
Optimization of DG in the system. A comparison between the Firefly Algorithm (FA) with 
Power losses Genetic Algorithm (GA) is also demonstrated in this paper. The results 
obtained have proven that the Firefly Algorithm has a better capability at 

improving both the voltage profile and the power losses in the system. 
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1, INTRODUCTION 

In a distribution network, the distributed generation are used for several purposes, for instance 
reducing power losses, obtaining better voltage profile along the feeders, as well as to maximize the 
increment of transmitted power in transformers and cables. For power loss reduction, the distributed 
generation which have been installed in power distribution systems are used for compensation of reactive 
power. However, the distributed generation (DG) installation in distribution systems requires consideration of 
their suitable placement and sizing for distributed generation. The main objective of this work is to find the 
optimal location and sizing of the distributed generation for in order to minimize the losses of real power and 
also for the improvement of voltage profile. However, the installation of distributed generation units at non- 
optimal locations and also non-optimal size can causes several problems such as voltage flicker, damaging 
the voltage state, system losses, harmonics and power system stability which it can lead to an increase in 
operational costs [1-2]. 

Therefore, the optimization method used in recent studies shows that it is able of showing the 
suitable answer for a particular distribution system can be very suitable in power system development. 
Studies shown that, various optimization method have been created for solving the ideal capacity and 
placement of distributed generation in the distribution networks. The optimization technique included such as 
genetic algorithm, tabu search, heuristic algorithms [3], analytical based methods [2] and metaheuristic 
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algorithms which were developed based on swarm intelligence in natural world such as Particle Swarm 
Optimization (PSO), Ant Colony Optimization as well as Artificial Bee Algorithm techniques [4]. There have 
been numerous work carried out by researchers worldwide in the area of locating the optimal capacity and 
placement of DG using all these optimization techniques. In the work done in [12], the optimal size and 
location of DG was carried out on an IEFE 69 bus test system and it was found out that the Firefly Algorithm 
was able to minimize the total power losses in the system. In another work, the optimal allocation of several 
DGs utilizing a combination of analytical and genetic algorithm methods to minimize the system losses was 
compared using the 33-bus and 69-bus IEEE test distribution systems. The simulation results show that the 
proposed method resulted in lower power losses obtained [13]. Other than that, in the work done in [14], a 
method known as gravitational search algorithm (GSA) was used to minimize the losses and total harmonic 
distortion in the system while locating the suitable placement and sizing of DG. The results have shown that 
GSA performs better than PSO in having better convergence rate and fitness function as well as minimizing 
losses in the system. On the other hand, the optimal location, sizing, and contract pricing of DG in electric 
distribution systems were found using three optimization methods of Scatter Search, GA and PSO. It was 
found that that GA and PSO performed better in solving the multiobjective problem [15]. 

In this research work, the voltage and system losses analyzation with the incorporation of distributed 
generation (DG) in the system is carried out using a metaheuristic optimization technique known as the 
Firefly Algorithm. This algorithm which is inspired by nature and it is a metaheuristic optimization technique 
which works based on the public (sparkling) nature of fireflies. Fireflies are considered by their flashing light 
created by natural activity known as bioluminescence. In this cases, the flashing of light may perform as the 
key of courtship signals for mating with other fireflies and also be used for giving a warning to potential 
predators. This flashing light will be connected with the objective function to be improved, which will create 
the potential to express new optimization algorithm [5]. 

This paper presents a simulation of load flow analysis of an IEEE 33-bus radial distribution network 
for finding the losses of power and voltage improvement. The data for IEEE 33-bus distribution network is 
taken from [11]. In this study, the firefly algorithm (FA) has been proposed in order to determine the best 
answer for the ideal allocation and size of distributed generation (DG) in an IEEE 33-bus distribution test 
network. The simulation of the test system was carried out in MATLAB software. The paper is structured as 
follows in which Section 2 presents the problem formulation for this work, the optimization algorithm and 
methodology for the firefly algorithm. Next, the results and discussion obtained from the simulation work is 
discussed in Section 3. Lastly, in Section 4, the conclusions from the research work is presented. 


2. PROBLEM FORMULATION 

In this paper, the main objective function of this optimization problem is aimed at obtaining loss 
minimization and improvement of voltage profile. As the losses is minimized, the voltage profile at every bus 
is estimated to be improved. The equations of the optimal location of distributed generation is stated as 
below: 


X = [Xqp, X15) Xnp Xnz] (1) 


where S shows the size of distributed generation while L shows the location of distributed generation and n is 
the number of distributed generation component that will be connected in the network. Equation (1) is 
included in the load flow data and then it 1s executed for obtaining the total losses in the system. The 
objective function, f(x) is the outcome of the sum of losses in the network as the Ploss is minimized. The 
following equation is formulated based on equation (1) which is the objective function for this work: 


f (2c) = min(X}¢ Ploss) (2) 


where line is referred to number of lines in the distribution network. Ploss can be obtained by using equation 
(3) and is also based on the parameters in equation (1). 


Ploss = Yihy Dyes laij( PP; + Q:Q;) — Bij(QiP; — P:Q;)) (3) 


For this research work, firefly algorithm has been introduced as an application for solving the ideal 
placement and capacity of distributed generation problems. As mentioned in the introduction, firefly 
algorithm is one of the swarm based optimization algorithm where it is performing based on the sparkling 
performance of all the fireflies. The steps or description of all parameters required in modelling the algorithm 
for finding the ideal size and location of distributed generation is presented below. Based on the founder of 
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firefly algorithm, he stated that [5]: 

1. A single firefly will be lured to other fireflies according to their sex even though there are unisex. 

2. The brightness is proportional to their attractiveness. Therefore the firefly which is less bright will move 
towards the brighter firefly. As there are proportional to each other, their attractiveness will decrease as 
their distance increases. Hence, fireflies will move randomly if there is no brighter firefly. 

3. The landscape of the objective function will determine the brightness of a firefly 

It can be assumed that the attractiveness of a firefly is observed with its brightness where it is related 

to the objective function. For the simple case of the maximum optimization problems, the attraction which is 
‘B” is virtual where it can be seen in the eyes of the other fireflies. Hence, the distance rj; will be exchanged 
with firefly ‘j’ and firefly ’1’ [5]. There are three important elements that have one has to follow in order to 
find the optimal placement and sizing of distributed generation in the network using firefly algorithm. The 
parameters were selected for simulation when using the firefly algorithm technique. These include the 
attractiveness By = 1.0, ae [0,1] and the attractive quality or the ingestion coefficient y=1.0, which are used 
to ensure a rapid interaction of firefly algorithm for the calculation to the ideal solution. 


2.1. Attractiveness 

In the firefly algorithm, the application of the formulation of attractiveness and light intensity is very 
important in this algorithm. Light strength is proportional to the attraction of all the fireflies and therefore, 
the attractiveness function is stated as follows: 


By = By * exp(—yj}),withm > 1 (4) 


where, Bg is the attractiveness at r = 0 while r is the separation between several fireflies. Furthermore, y is an 
ingestion of constant which dominates the diminishing of light power. 


2.2. Distance 
A Cartesian equation is taken after deciding the separation between any two fireflies 1 and j, at 
positions X; and X;: 


= is 7 x;l| = ADyalesr: — XK) %2 (5) 


where d is the quantity of measurements, for d = 2, and xj, is the kth part of the spatial arrange x; of the ith 
firefly. Therefore, rj; which is given as: 


ry = (xi — x) + (yi - yi) (6) 


Moreover, the computation of the separation of r can likewise be characterized by utilizing other 
separation measurements which depends on the idea of the issue for instance Mahalanobis separation or 
Manhattan distance. 


2.3. Movement 


The firefly movement j which is interested by a more attractive or the brighter firefly 1 is given by 
the stated formulation below: 


X; = Xi + Bo * exp(—yj;) * (x; — X; ) +a * (rand — =) (7) 


A procedure for implementing the firefly algorithm with the load flow analysis of IEEE 33-bus 
radial distribution network in order to find the optimal location and sizing of distributed generation in the 
system are summarized in Figure |. After all the solutions are obtained, the results are discussed in part 3. 


Optimal sizing and location of distributed generation for loss minimization ... (Muhamad Najib Kamarudin) 


424 O ISSN: 2502-4752 


Start Firefly Algorithm 
ae 


Insert the firefly param eter, line and 
bus data. min and max size of DG 


= 

Run power flow and determ ine the 

power losses and voltage at each bus 
= 
Com pute the relevant form ula for determ ining 

tge optimum DG sizes at each bus 

= 
Calculate candidate of DG 
size and represent as fireflies 

= 


Print the size and location of 
Initialize firefly offspring i DG 
= a, 


Yes 
Generate 
random firefly j Iteration 
= maximum 
Sum and com pare their 
baal 


No 


fitness (minim um loss) 


Check if1<j Determine the best fitness 
=at f 





Figure 1. Flowchart for finding the of ideal sizing and location of DG using firefly algorithm 


3. RESULTS AND ANALYSIS 

The six firefly algorithm parameters variable values used are given below in Table 1. The simulation 
study was carried out by determining the system loading condition. The IEEE 33-bus system selected for 
base kV and MVA is 12.66kV and IOOMVA. The data for the line and load for the IEEE 33-bus system was 
taken from MATPOWER 6.0. The test system incorporated with the optimization algorithm had been 
simulated on MATLAB R2018a software and run on MACBOOK with IOS operation system. Figure 2 
below shows the 33-bus radial distribution network. The sum of the load of the network is 3.72 MW for the 
active power and 2.3 MVAR for the reactive power while the active power loss 1s 281.58 kW. 


Table 1. The Firefly Algorithm Properties Used In the Optimization Formulation 


Firefly Algorithm Parameters Values 
Population of Firefly 20 
Highest number of iteration 30 
Number of DG 1 and 2 
Capacity of DG 0.01MW to SMW 
Alpha 0.25 
Beta 0.2 
23 24 25 


26 27 28 29 30 31 32 33 
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Figure 2. The IEEE 33-Bus radial distribution network 
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To further investigate the effectiveness of the firefly optimization algorithm in finding the ideal 
sizing and location of the distributed generations, the following three cases are simulated: 


Case 1: Simulation of the system without DG, to determine the base case scenario. 
Case 2: Simulation of the system with one DG included. 
Case 3: Simulation of the system with two DGs included. 


In order to show the advantages of using firefly algorithm, the three simulation cases are also carried 
out using another optimization method which is Genetic Algorithm. Table 2 below shows the outcomes of 
the simulation of the power flow of IEEE 33-bus distribution network using the MATLAB software 
integrated with the optimization algorithm for finding the optimal sizing and placement of the distributed 
generation using the two optimization algorithms. 


Table 2. Comparisons of Results between Firefly Algorithms with Genetic Algorithm for the Three Case 


Studies 
Total Real Power Hus vouase NiicraRR bus Optimal size 
losses and location 
(MW) Bus number Voltage (p.u) Size (MW) Location 

Case 1(without DG) 0.193 18 0.8820 - - 
Case 2 (with one DG) FA 0.178 18 0.9189 1.5214 7 
GA 0.219 18 0.9354 1.4516 7 
Case 3 (with 2 DGs) FA 0.134 18 0.9571 0.5771 7 
2.0871 27 
GA 0.202 18 0.9619 0.5645 7 
2.0968 27 


Based on the results obtained from Case 2 which is by connecting one distributed generation into the 
system, the network losses recorded is 0.178 MW and the minimum bus voltage obtained is 0.9034 V. This 
power loss result is better than Case | which is the base case where in the base case, the total real power loss 
obtained was 0.193MW. The simulation also shows that after installing a single distributed generation into 
the system, the optimal location recorded is at bus number 7 and the optimal sizing of the DG is 1.5214 MW. 
Case 3 considered connecting two distributed generations into the network and the network losses recorded is 
0.134 MW and the minimum bus voltage obtained is 0.9371V. It can also be seen that the suitable locations 
after installing two distributed generations are at bus number 7 and 27 with the appropriate size of 0.5645 
MW and 2.0968 MW. It can be seen from the results obtained using Genetic Algorithm (GA), the ideal 
placement of distributed generation is similar to the ones obtained using firefly algorithm. For case 2, the 
optimal location obtained using GA for installing the DG is obtained to be at bus 7 and for case 3 the ideal 
placements obtained are at buses 7 and 27 respectively. 

The total real power losses recorded by firefly algorithm is better than achieved by using genetic 
algorithm for all the three case studies. In case 2, by installing one DG in the system, it can be seen that the 
total real power loss is reduced by 7.8% compared to the base case without distributed generation installation. 
Next, for case 3, reduction on the total real power losses obtained is 30.6% compared to the base case. The 
calculations for obtaining the percentage of reduction of power loss is shown next: 

(1 DG) — Real power losses = (0.193-0.178)/0.193 x 100% = 7.8% 
(2 DGs) — Real power losses = (0.193-0.134)/0.193 x 100% = 30.6% 

Figures 3 and 4 below shows the voltage profile obtained with and without the inclusion of DG in 
the system. It is seen in Figure 4 that with the installation of DGs, the voltage profile is better and the curve is 
within the optimal value of 1.0 p.u. 
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Figure 3. Voltage profile without DG included in Figure 4. Voltage profile with one and two DGs 
the system included in the system 


4. CONCLUSION 

This paper has presented the implementation of firefly optimization algorithm in finding the most 
suitable placement and capacity of DG when DG 1s incorporated in the system. It can be seen from the results 
obtained, the firefly algorithm is better compared to another well known optimization method of genetic 
algorithm (GA) for determining the optimal placement and capacity of distributed generations in the 
distribution system. The results show that integrating the distributed generation with the optimal capacity and 
the suitable placement in the distribution network is able to reduce the stil Sf line power losses and the 
improvement of voltage profile. It is also seen that the total power losses and the total voltage deviation of 
the network is better with the integration of two DGs compared to the case with one DG. It is suggested that 
for future work, an improvement in the problem formulation and objective function could be achieved. In 
conclusion, the metaheuristic Firefly Algorithm (FA) has shown to be promising in determining the optimal 
placement and capacity of the distributed generation while a the total real power loss and 
improving the voltage profile. DG 
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